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Abstract 

The different natures of approximate symmetries and their corresponding first inte- 
grals/invariants are delineated in the contexts of both Lie symmetries of ordinary 
differential equations and Noether symmetries of the Action Integral. Particular note 
is taken of the effect of taking higher orders of the perturbation parameter. Ap- 
proximate symmetries of approximate first integrals/invariants and the problems of 
calculating them using the Lie method are considered. 

1 Introduction 

According to Gazizov |^] the first publication devoted to approximate symmetries was by 
Baikov et. al. Q (see also |^]). The concept of an approximate symmetry goes back to 
at least the late seventies [|l2| when the observation was made that, since the differential 
equations which arise in mathematical modelling are invariably approximate, one should 
in fact be considering approximate symmetries. The observation may have been a conse- 
quence of a number of papers on applications in Cosmology in the sixties and seventies 



(cf. 14, 21, 26, These applications continue pq]. More recently in the context 



of differential equations the concept was regarded as sufficiently important for Head |1C| 
to add a procedure for the calculation of approximate symmetries to his Program LIE 
in 1997. The procedure is for computations to the first order in the perturbation param- 
eter e, but the program is easily adapted for the computation to be made to any order 
(subject to the capacity of the machine). 

We recall that a differential equation has a Lie symmetry if there exists an infinites- 
imal transformation which leaves the differential equation unchanged. The infinitesimal 
transformation 

y = y + ar], 
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where a is the infinitesimal parameter, x is the independent variable and y the dependent 
variable (the extension to more than one of each is simple enough in principle), can be 
written as 

x = (l + aG)x, y = (l + aG)2/, (1.2) 

where G is the differential operator 

G = + ^idy. (1.3) 

To deal with differential equations the operator G must be extended to give the infinites- 
imal transformations in the derivatives induced by the infinitesimal transformations in 
the independent and dependent variables. The nth extension of G, which deals with all 
derivatives up to the nth derivative, is [20| 

n 

+ (1.4) 



where 



dx dx 



(1.5) 



An nth order differential equation, 

i?(x,y,y',...,yW) =0, (1.6) 
is invariant under the infinitesimal transformation generated by G if 

= 0. (1.7) 

|£;=o ^ ' 

Each differential operator G, with this property is called a symmetry of the differential 
equation and the aggregation of all such symmetries forms the Lie algebra of the dif- 
ferential equation. A first integral/invariant, I{x,y,y' , . . . ,y^^~^^), is associated with a 
symmetry G if 

= and ^ = 0. (1.8) 

dx 

Instead of a differential equation one can consider the Action Integral as the object 
to maintain its invariance under an infinitesimal transformation. In this case we have 
Noether's Theorem ||2^ which states that G is a Noether symmetry of a Lagrangian, 
L (x, y,y' , . . . , y^"^) , if there exists a function /, such that 

/' = G["1l + ^'L. (1.9) 

The existence of the symmetry is independent of the Euler-Lagrange equation. When 
the Variational Principle is invoked to give the Euler-Lagrange equations, (|1.9| ) can be 
manipulated to give the first integral/invariant 

I = f-[^L+{v,-y',()]^. (1.10) 
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All of these considerations have been made without reference to the nature of the coeffi- 
cient functions ^ and ij. In Lie's original work on the symmetries of differential equations 
he considered both point [17| and contact transformations p^] while Noether used gener- 



alized transformations ab initio. There is in fact a priori no restriction on the functional 
dependence of the coefficient functions and increasingly use is made of generalized and 
nonlocal symmetries to expand the classes of problems which can be solved. If algebraic 
properties are the primary consideration, it is better to stay with point or contact symme- 
tries, otherwise there is no need for constraint. In the end the choice of type of infinitesimal 
transformation to be used is determined by a utilitarianism appropriate to the occasion 

No matter the type of infinitesimal transformation under consideration, it is always 
possible to introduce the approximate form of such a transformation if the differential 
equation under investigation contains a small constant so that the equation can be re- 
garded as a perturbation of the equation with the constant set at zero. Suppose that the 
differential equation is 

E(x,y,y'...,y^''\€)=0, (1.11) 



x,y,y ,...,y' 

where e is the small constant and should not be confused with the a used to denote the 
infinitesimal parameter in (|l.l| ) and (|l]^) . The generator of an approximate infinitesimal 
transformation of order k can be written as 

k 

G = ^e^Gi, (1.12) 

i=0 

where 

-^ = 0. (1.13) 
The generator ( 1.12| ) will be an approximate symmetry of ( |l.ll ) if 



= O (e'^+ij . (1.14) 



We note that, when n = in ( 1.11 ), we have the case of a function and the definition of 
an invariant associated with a symmetry group is generalized to that of an approximate 
invariant or, more correctly, an invariant function of an approximate symmetry group. 
(In the language of jet bundles, cf. Cotsakis, Leach & Pantazi one is extending the 
"approximate" vector field from one on the manifold — in the case of a function — to that 
on the re*'* jet bundle of the total space. A differential equation is then a variety living in 
the n^^ jet bundle and in our case we may reinterpret this as the unperturbed entity of 
the e-perturbed family of varieties related to the approximate symmetries.) 

Noether's Theorem can be adjusted as follows. We consider transformations generated 



by symmetries of the form ( 1.121) so that (|1.2| ) becomes, to order k in e. 



x + a^e^Ci, y = y + a'^e'r]i, (1.15) 



where 

Gi = Cidx + i]idy, 



(1.16) 
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From Noether's theorem we have that (|1.12|) is an approximate Noether symmetry of the 
action integral 

A= L{x,y,y',a)dx (1.17) 

if it leaves the action integral ( |1.17 ) unchanged, i.e. 

/ L{x,y,y,a)dx= L {x,y,y' ,a) dx + aF + O (e'^+^V (1.18) 
Jxo Jxo 

where aF is the infinitesimal contribution from the boundary terms and, naturally, is 
written in terms of an expansion in e as 

aF = aj:^F., § = = §. (1.19) 

i=0 

The use of ( |1.15| ) in ( |l.lg| ) up to order leads to the Killing-type equation [^] 

k k 



1=0 i=o ^ y J 

+w t+ 6^.} +0 



(1.20) 



where prime denotes total differentiation with respect to x, we have written L = Lq + eLi 
in an obvious notation and all terms of order greater than are neglected. Hence ( 1.1(]| ) is 
an approximate Noether symmetry of ( 1.17 ) if the Killing type equation ( |l.2[)| ) is satisfied. 

Nothing has been said about the dependence of the functions i]i and fi apart from 
their lack of dependence upon e. If one seeks approximate point transformations, these 
functions depend on x and y only. Otherwise the functions can be allowed to depend 
on X, y and derivatives of y. Indeed there is no reason why one should not contemplate 
nonlocal approximate symmetries. A theory of approximate Noether symmetries based 
on the concept of approximate infinitesimal transformations can be found in However, 
that theory is limited to the first order in e only and does not lead to the results presented 
in this paper. 

The rest of this paper, which is the third in a series |15, 16| devoted to the investigation 
of the connections between the three main topics in dynamics, viz symmetry, singularities 
and integrability, topics which superficially are unrelated in their mathematics and yet 
are intimately intertwined, is organized as follows. In the next section we examine three 
simple second order ordinary differential equations for approximate symmetries. The 
equations are 

(i) the simple harmonic oscillator treated as a perturbation of the free particle, 

(ii) the Ermakov-Pinney equation treated as a perturbation of the equation for the 
simple harmonic oscillator and 

(iii) an autonomous Emden-Fowler equation of order two considered as a perturbation 
of the free particle equation. 
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Each of these equations illustrates certain aspects of the properties of approximate 
symmetries. We have chosen simple equations so that the properties of the symmetries 
will not be obscured by complexity of the equations. In § 3 we consider the properties of 
approximate Noether symmetries applied to the Lagrangians of the differential equations 
above. In § 4 we examine some of the problems which arise in the determination of ap- 
proximate first integrals/invariants using the Lie method and the approximate symmetries 
associated with them and in § 5 we present our observations and conclusions. 



2 Approximate Lie symmetries: Three case studies 

2.1 The simple harmonic oscihator 

We write the differential equation as 



X + ex = 0, 



(2.1) 



where e is a small parameter. Naturally this is somewhat artificial since any value can 
be used as the coefficient of x by the simple expedient of rescaling time. We recall that 
the number of exact symmetries is eight no matter the value of e. The calculation of the 
symmetries is straightforward and we list the approximate symmetries at 0(e), O (e^) and 
O (e^) . They are 



0(e): 
Gi = xdx, 



Gs = ( 1 - -et' ] d,, 



G2 = e{2tdt + xd^), 
1 



Ga 



G5 = x[t- ^etA dt + x^(l- l-eA d„ Gq = (t^ - ^eA dt + x ( t - ^th ] d„ 



3! 



G7 = x{l--et^]dt-txhd,, 



2! 



Go 



'l-2et^)dt-2txedr,, 



Gil = edt, 

Gi3 = etdt, 

Gi5 = e {txdt + x'^dx) 



Gs = exdt 



Gio= (^t-^et^^dt + x{l-t^)d,, 

Gi2 = edx, 

Gi4 = etdx, 

Gi6 = e {t'^dt + txdx) ; 



(2.2) 



Gi = xdx, 



G, 



1 - l^t' + ^eh^ ) dx, 



G2 = e2 {2tdt + xdx) , 

1 o 1 

3!' +5! 



G, 



t 
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G. = e 



Ge=[t'- let' + ^eh'^ dt + x^t- ^et' + ^^e't' ] d., 
Gr = x[l- \et' + j/A d, - x' (et - ^^ehA 



Gg = ( 1 - 2et^ + ^e^t' ] dt - x { 2€t - ^e^t^ ] d^, 



Gio — e 



X ( 1 — l-et^ ] dt — etx'^dx 



t - ^et' )dt + x{l- t^) d. 



1 



Gil = e^dt, 



Gi2 = e 1 - -et' d, 



Gi3 =\et- -eV I dt - eH'xd, 



GiA = ^\^- ^et^j dx, 
Gig = {t^dt + txd^) , 

Gu = e (^t' - ^et') dt + ex (t' - ^et^ ) d., 



Gi5 = (txdt + x'^dx) , 



Gis = x{t- ^et^ + ^eH' ] dt + x^ 1 - '-et^ + ±,eh' ) 9., 



Gi9 — e^dx, 

G21 = ehdt, 



G23 = e[{l-2et^)dt-2etxdx] , 



1 2 1 
2^' + 1- 

G20 = e^xdt + exdx, 

G22 = ^'^tdx, 

G24 = e [(1 - 2et^) dt - 2etxdx 



(2.3) 



O (.3) : 
Gi = xdx, 

G3-- 



1 _ let2 + 1,2^4 _ 1 3^6 ) Q 
2 4! 6! ' 



G2 = e^ {2tdt + xdx) , 



Ga 



G5 = e' 



1 



3! 



x\t- -er ]dt + x' 11- -et' d, 



1 



2! 



1 



t - ^et^ + ^e'^t^ - ^e^t' ) d.. 



.3.7 



3! 5! 



7! 



1 



45 
1 



GQ = e[t' - -et"" + — e^t^ ]dt + ex[t- -ef' + — e^t^ a 



15 



G7 = X ( 1 - ^ef + i-eH' - -eH^ \ dt - x' \ et - -e'f + -eH'> \ d, 



.3+6 



G8 = e' 



6! 



1 - ^et^ ] dt - x'^etdr 



21 4! 
1 
2 



.2+3 



.3^.5 



3! 



5! 
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Crg — 


( 1 o j-2 , 2 2,4 

1 1 - 2et + 3^ * " 


Gio = 


/ o 1/1 2 o 

r-r^ +45^ 


Gii = 




Gi3 = 




Gi5 = 


[txdt + x^d^) , 


Gn = 


/ 1 o 1 

xt -xet'^ + —: 

V 3! 5! 



^-^.3,s),^^(,._l,.,,.^l,.,.,4_l, 

G16 = (i^^i + txd^) , 

- ^xet-^ + ^xe^t^ - ^xe^t'^ 1 dt + ( xt - -xet^ + —xe'^t^ —xe^f ] d^, 

,3! 5! 7! y V 3 15 315 ' ' 

Gig = e^d,, G20 = (1 - ^e*') 

G21 = ehdu G22 = ehd,, 

G23 = (e^ - 2e^t^) dt - 2x€^td^, G24 = (^2e^t - ^e^t^) ^ (^^^2 _ 2xt^e^) 

G25 = (2et - ^eh^ + ^eh'] dt - ( 2xeh^ - Ixeh^ \ 

G2, = {t- \et^ + - ^eh-^^ dt + [^^xeH' - xet" - ^^xeh^ ) 

2.2 , 2 3 4\ „ / 0^,2. 4 , 3 



G27 = e - 2tH^ + -e'^t* ]dt-[ 2xeH - -xe^^ 9^, 



G28 = [xe - -xeH^ + -^xeh^j dt - i^xhh - ^xhh^ ] 9,, 
G29 = (^eh' - ^eh'^ dt + (^xeh - ^xeh'] d. 
G3o=(^et-^eh^ + yh'^d^, 

G31 = (xrf - ^xeh^ + ixe^ts) + (x^e - ^x^e^t^ + ±^2^3^4 , 
G32 = e'xdt. 



(2.4) 



We observe that there are sixteen approximate symmetries at 0(e), twenty-four at 
O (e^) and thirty-two at O (e^) . On closer inspection we see that there are different cat- 
egories of approximate symmetry which would not become apparent if we worked only 
to 0(e). The homogeneity symmetry, xd/dx, which is common to both perturbed and 
unperturbed equations occurs at all orders of e. Other symmetries occur at levels of in- 
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creasing approximation, for example G3, in which one sees the role of the pertm'bation 
parameter bringing additional terms of a MacLaurin expansion. We recognise the initial 
terms of the expected sine and cosine series. Then also we see, again with the example of 
G3 at 0(e), that G12 is simply e times the first order G3. Hence G3 is a true approximate 
symmetry in contrast with, say, G12 which is simply a truncation symmetry. At 0(e) one 
cannot distinguish between the natures of the two types of symmetry. Consequently we 
can recognise three types of "approximate" symmetry. The first category comprises those 
symmetries which are common to both perturbed and unperturbed equations and which 
occur at all orders of the perturbation parameter. The second class are those symmetries 
which are truely approximate as one can see the increasing level of approximation with 
increasing powers of e. They also occur at all levels of approximation. The final category 
are those which exist by virtue of the power at which the approximation was truncated. 
They can also have expansions of a form similar to the second category, but are distin- 
guished by never having a term which is independent of e. The number of approximate 
symmetries increases by eight with each additional order of approximation. 



2.2 The Ermakov— Pinney equation 

The Ermakov-Pinney equation is one of those equations which appears in divers 

areas and is intimately associated with the sl{2,R) subalgebra which is the characteristic 
of all scalar nth order ordinary differential equations of maximal symmetry |20|. We write 
the equation as 



x + uj'^x = At (2.5) 



which has the implication of a perturbation from the equation for a simple harmonic 
oscillator. We note that both uj and e can be scaled to any numbers by rescaling of the 
variables x and t. The Ermakov-Pinney equation has three exact symmetries so that here 
we are looking at the effect of a decrease in exact symmetries from eight to three. Again 
we list the approximate symmetries at 0(e), O (e^) and O (e^)- They are listed in (|2.6| ). 

The pattern in ( ^.61 ) is perhaps easier to discern than in (p.2|), (|2.3| ) and (|2.4| ). The first 
three approximate symmetries, Gi, G2 and 6*3, are repeated at 0(e), O (e^) and O (e^) as 
the order of the approximation is increased and account for the first twelve approximate 
symmetries (where appropriate). These symmetries, with the algebra sl{2,R), are the 
symmetries common to both the equation for the simple harmonic oscillator and that 
for the Ermakov-Pinney equation. The remaining symmetries are all at the maximum 
power of e applicable and so are a consequence of the truncation of the perturbation 
expansion. They represent the remnants of the homogeneity symmetry (G15 and the 
solution symmetries (Gie and G17) of the original simple harmonic oscillator. The number 
of approximate symmetries increases by three at each additional order of approximation. 
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0(e) O (e2) O (e3) 

Gi =dt Gi = dt Gi = dt 

G2 = sin 2LLitdt G2 = sin 2ujtdt G2 = sin 2ujtdt + lux cos 2ujtdx 
+ajx cos 2ujtdx +ujx cos 2ujtdx 

G3 = cos 2ci;t9t G3 = cos 2ujtdt — (jJX sin 2u)tdx G3 = cos 2Li;t9t — lux sin 2u!tdx 
—ujx sin 2ujtdx 

G4 = e9t G4 = e9t G4 = e9t 

G5 = e sin 2a;<9t G^ — e [sin 2wt9f + wx cos 2a;<92:] G5 = e [sin 2wt9t + wx cos 2ujtdx\ 
+eojx cos 2ujtdx 

Gq = e cos 2a;t9t Gg = e [cos 2(^^94 — sin 2a;<9a;] Gg = e [cos 2ujtdt — ujx sin 2a;<92;] 
—eujx sin 2(jjtdx 

G- = e^dt G7 = e^a, 

Gs = sin 2ojtdt + e^ujx cos 2ujtdx Gg [sin 2wi9t + lux cos 2ujtdx 

Gg = [cos 2ujtdt — lux sin 2wt9a;] Gg = [cos 2ci;t9t — wx sin 2ujtdx 



Gin = e^a, 



Gil = e'^ [sin 2ci;t9t + ujx cos 2cjt9:i.] 
G12 = e'^ [cos 2ujtdt — wx sin 2<x't92;] 
Gi3 = exdx Gi3 = 6^x9^ G13 = e^xdx 

Gi4 = esincijtc^j. G14 — sinujtdx G14 = sintotdx 

Gi5 = e cos ujtdx Gi5 = e^coswtc}^; G15 = e'^ cos wtc^a; 

G16 = exsinwtSt Gig = [xsina;t9t + wx^ coswtc}^.] Gig = e'^ [xsinwtSt + ljx'^ cos uutdx] 
+eujx'^ cos wtSa; 

G17 = ex cos ujtdt G17 = [xcosw^Sj — wx'^ sinwt^a;] G17 = e^ [xcoswt^t — uox"^ sinwt^a;] 
— ewx^ svaujtdx 



(2.6) 



2.3 Emden— Fowler equation 

We consider the Emden-Fowler equation 

X = ex^ (2.7) 
i.e. as a perturbation of the free particle equation. Again the parameter can be set at any 



value by rescaling. We find the approximate symmetries as are listed in (2.9). 

The symmetries Gi and G2 are repeated at higher orders of approximation and so are 
common to both equations. In the normal listing of the Lie point symmetries of the free 
particle one would have the two symmetries 

Xi = xd^ and X2 = tdt + ^xd^, (2.8) 

the homogeneity symmetry and the middle element of the sl{2, R) triplet in the form ap- 
propriate to a second order ordinary differential equation [^] , but any linear combination 
of Xi and X2 is also a symmetry. The combination which gives G2 is common to both 
equations and so is maintained at all orders of e whereas the combinations which give Gg 
and G12 are not inherited, but are simply a consequence of the order of the approximate 
symmetry. We also observe that the perturbations of Gi occur to the order of e being 
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used in the approximate symmetry whereas those of G2 are at one order less. This is be- 
cause the differential equation, (^), is used in the calculation of the second approximate 
symmetry but not the first which simply annihilates the equation. 



0(e) 
Gi = 
G2 = 

G3 = 



dt 
tdt 



2xdx 



Gs — 

G9 = 

Gio = 
Gil = 
G12 = 

Gi3 = 
Gl4 = 
Gi5 = 

+ {6t 
G16 = 
+3(1 



etdt 
£xdt 
: etdx 

- £xdx 

- e 



txdt 



txdr 



f 2et^x) dx 
+ et^x) dx 



Gi =5* 

G2 - tdt - 
G3 = edt 
G4 = e [tdt 
G5 = e^dt 



Gs = e^dx 
Gg = eHdt 



2xdx 
— 2xdx 



e xdt 
eHdx 
e^xdx 
txdt 

edt^ 

Gi5 = e {^t^dt - 



Gio 
Gil 
G12 
Gi3 = 
Gi4 = 



- x^dx 
txdx 

{6t + 2et^x) dx 



2xdx 
2xdx] 



Gi^dt 
G2 - tdt 

G3 = edt 
Gi = e [tdt 
G5 = e^dt 

Ge = [tdt - 2xdx] 

Gj = e^dt 

Gs = e^dx 

Gg = eHdt 

Gio = e^xdt 

Gn = ehdx 

G12 = e^xdx 



Gi3 ~ 
Gi4 = 
G15 = 



txdt 

t'dt 



- x^dx 
txdx 
(6H 



2et^x) dx\ 

G16 = e \^tt^dt + 3(1 + tt^x) dx\ G16 =- <? \2d?dt + 3 (l + d?-x) dx\ ■ 

(2.9) 



3 The Noether symmetries and approximate 
first integrals /invariants 

3.1 The simple harmonic oscillator 

The simple harmonic oscillator has the Lagrangian 

L = \{±'- ex'] . (3.1) 

Given the observations made in the previous section we know the pattern of the behaviour 
of the symmetries and so we will simply list the results for the first two orders of e. We 
consider point symmetries only. Because of the direct relationship between symmetry and 
first integral we list the Gs, fs and Is together. The results are given in (3^) and ( |3.3| ). 



The only Noether symmetry of the free particle which is inherited by the simple har- 
monic oscillator is Gi. Since both Lagrangians have five Noether point symmetries, at 
0(e) we see that there are two symmetries due to the inheritance of the exact symmetry 
Gi, four approximate symmetries of the free particle which persist as truncation symme- 
tries and four genuine approximation symmetries of the simple harmonic oscillator. We 
see the beginnings of the MacLaurin expansions for the sine and cosine functions. The 
behaviour of the integrals reflects these distinctions. At O (e^) we gain an extra five sym- 
metries. One of these is due to the inherited symmetry and the other four because the 
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four approximate symmetries of the simple harmonic oscillator at order one will have the 
same behaviour as the inherited symmetry at order two. The distinction between the two 
classes is that for the inherited symmetry all occurrences are simply multiples of the ze- 
roth order symmetry by powers of e whereas the basic approximation symmetries increase 
the degree of approximation with increasing powers of e and the others are multiple of 
approximate symmetries at lower order. The truncation symmetries do not multiply, but 
simply persist at the order of the approximation being used. 
To the first order in e we have 



Symmetry 


Gauf 


function 


Invariant 






/i = 





/i = i [x' + ex') 




G2 = edt 


./2 = 





h = \ex^ 

I3 = X {tx — x) — et (x^ - 




G3 = [2t - ie(2i)3] dt 
+x [2t - ^e(2i)2] d. 


/3 = 


-2etx'^ 


■ 2txx + It'x'' 










Gi = ei2tdt+xd,) 
G5 = (1 - ^et^) d. 


.h = 





/4 = ex {tx — x) 




/5 = 


—etx 


4 = (1 - ^et')x + etx 




Gq = edx 


/6 = 





h = ex 




Gj^{t- l,et^) d. 


/7 = 


(1 - l,et^) X 


h^{l-^ee)x^{t- 


^et^) X 


Gs = etd^ 


./8 = 


ex 


Ig ~ e {x — tx) 




G,^[l,{2tf-^e{2tf]dt 
+x [2t - ^e(2t)3] d^ 




(l-i6(2t)2)x2 


Iq = {x- txf - let' (3a; 


— tx) {x — tx) 




/lo = e{x — txf 




Gio = ^e{{2t)^dt+x{2t)d,) 


./io = 







(3.2) 



and to the second order in e we have 



Symmetry 


Gauge function 


hivariant 


Gi =9* 


/i=0 


/i = 1 (i"* + ex"*) 


G2 = edt 


/2=0 


h = \e {x^ + ex^) 
h = lex^ 


G3 = e'dt 


/3 = 


G^^[2t-^e{2t)^ + ^e^{2t)^]dt 
+x [2t-^e{2tf + y\2t)^]d, 


U = -ex^ [2t - ^e{2t)^] 


I4 = X (tx — x) 


-iet (3i2 -6tei + 2*2x2) 
+ ^6^*2 (5a;2 - 5txx + 3t2i2) 






Gs =e[2t~ iTe(2i)3] dt 
+ex [2t - ^e(2t)2] d^ 


/s = -2eHx^ 


/s — ex (tx — a;) 




-eH (x^ - 2tei + fi^i^) 


Ge = e^ {2tdt + xd^) 

Gy = (1 - irf^ + ^eH^) a. 


/6 = 


/e = e^i (ii — x) 


/7 = -ex {t - iet3) 


+ex [t - ^et^) 


Gs - e (1 - iet^) d. 


/s = -eHx 


h = e [(1 - ^,et^)i + etx] 


Gg = e^a^ 




h = e^x 


Gio = {t- ^et^ + ^£2^5) 


ho = {l-ht^ + ieH^)^ 


ho = (1 - ief2 + l^eH^) X 


-{t- i,6t3 + i,e2^5) i 


Gn =e{t- ^et^) d. 


hi = e (1 - k^t") X 


111 = e (1 - ^^et') X 
-e (t - -iei^) i 


G12 = eHd, 


fi2 ^e'x 


/12 (2; _ ij.) 


Gi3- [i(2t)2-ife(2i)4 + i,e2^6]g^ 


/i3 = (1 - M2t)^Y 


^13 = {x- tx)^ - 


+x [2t - y{2tf + ^eH^] d. 


+ ieH2t)^x' 


iet^ (3x - tx) (x - tx) 

+ j^e't'^ (5x2 _ 4^^^ _^ 2i2a;2) 


G,^^e{^,{2t)'~^e{2t)^)dt 
+ex {2t - ie(2i)3) 4 


/i4 = e (1 - hi^tf) 


-^14 = £ (x — tx)^ 




-ie2^2 „ - tx) 
/is = (3; — f^:)^ 


Gis = ^e^{{2t)^dt+x{2t)d,) 


/i5 = 



(3.3) 
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3.2 Ermakov Pinney system 

The Lagrangian for the Ermakov-Pinney system is 
• 2 ■? ^ 



L = l 

2 



2 2 
■ UJ X 



(3.4) 



To the first order in e we have 



Symmetry 



Gauge function 



Invariant 



G2 = edt 



/i-O 
.f2 = 



/i = i [x'^ + u^x^ + — 
h = |e (±2 +tj2^2^ 



G3 = sin 2(jjt9t + Lox cos 2ujtdx fs = —oj^x^ sin 2a;t = i ( — o^^x^ H — ) sin 2u)t 



G4 = cos 2ujtdt — Lox sin 2ujtdx 

G5 = e [sin 2Lutdt + lox cos 2ujtdx 

Ge = e [cos 2ujtdt — lox sin 2totdx 

G7 = e sin wiSj; 
Gg = e cos wtS^: 



/4 = a;^.T^ cos 2tji 

/s = — ecj^x^ sin2ijjt 

/e = ecj^x^ cos 2wi 

fj = ecja; cos 

/g = — eajx sintjt 



— tjxi cos 2iot 

Ii = \ (i^ — L?x^ + -^^ cos2tjt 
+ti;xi sin 2u:t 

Ir-, = ei (i^ — w^a;^) sin2[jt 

— eti;a;i cos 2Lot 

Iq = {x^ — uj'^x'^) cos 2(jjt 

+eijjxx sin 2[ji 

It = e{x sin — wa; cos wi) 

Is ~ e{x cos + sin ) 



(3.5) 



and to the second order in e 



Symmetry 



Gauge function 



Invariant 



/i =0 

/3 = 

/4 = — tj^a;^ sin 2ajt 
/s = w^x^ cos 2a;t 
Gg = e [sin 2ti;t9t + cos 2ujtdx\ /e = — ew^a;^ sin 2iot 



Gi=dt 

G2 = edt 
G3 = e^St 

G4 = sin 2ujtdt + wa; cos 2iotdx 
G5 = cos 2a;t9t — sin 2u)tdx 



h = \ [x^ + Lo'^x^ + — 
h = he fa;2 + w^a;2 + -1 



h ie^ (x^ +L0^x^ 



Lo'^X^ + 



sin 2Lot 



G7 = e [cos 2Lotdt — LOX sin 201^92;] 

Gs = sin ujtdx 
Gg = cos LOtdx 

Gio = e^ [sin 2ujtdt + ujx cos 2uotdx 
Gil = e^ [cos2ijjt(9f — U!xsm2u!tdx 



/r = ew^x^ cos 2u!t 

/g = e^tja; cos Lot 
/g = — e^tja; sintjt 
/lo = —e'^io^x^ sin 2a;< 

/ii = euj'^x'^ cos2iot 



■ uP'x^ + 



3 

/4 = I X 

— tjxi cos 2u}t 

h = \ {^^ - Lo'^x^ + cos2wi 
+Loxi sin 2tot 

Iq = if — uo'^x''^ + 1 sin2(jjt 
— ewxx cos 2a;t 
/. = (x^ 
+ecijx.x sin 2Lot 

/g = (.X sin Lot — wx cos tot) 
/g = (x cos Lot + LOX sin Lot) 

— e^wxx cos 2aj< 
III = (i^ - Lo'^x'^) cos2tjt 
+e^ti;xx sin 2wi. 



cos 2L0t 



\e^ fx^ - w^x^) sin2a;t 



(3.6) 
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We observe that the three symmetries, Gi, G2 and G3, of (|2.6| ) are common Lie and 
Noether symmetries of the equation for the oscillator and the Ermakov-Pinney equation. 
They constitute the algebra sl{2,R) which is the symmetry algebra for the latter and a 
subalgebra of the sl{3, R) algebra of the former. The solution symmetries are preserved by 
the order of truncation. The other three approximate Lie symmetries are not approximate 
Noether symmetries. The only serious first integrals/invariants correspond to these three 
symmetries. They are exact and not approximate. 



3.3 Emden— Fowler Lagrangian 

The Lagrangian for the Emden-Fowler equation (pT 



IS 



L 



1 



1 



-X + -ex . 



(3.7) 



To the first order in e we have 



Symmetry 



Gauge function Invariant 



Gi 


= dt 




h 


= 


h 


1-2 




G2 






h 


= 


h 




G3 


= 




h 


= 


h 


= ex 




G4 


= etdx 




h 


= 


h 


~ e{x — it) 




G5 


= e [2tdt -f 


- xdx] 


h 


= 


h 


= ex{x — it) 
= ^e{x - ti)"^ 




Ge 


= e [eOt 4 


- txdx\ 


h 


= 5^^' 


h 




G7 


= h^t^dt 


-f + \et^x) dx 


h 






= tx~e{^t^ 




Gg 


= let^dt + (1 + et^x) dx 


k 


= etx"^ 


h 


= X ~ e {tx"^ — 


t^xi + \t^i'^) 



(3.8) 



and to the second order in e we have 



Symmetry 



Gauge function Invariant 



Gi 


= dt 






h 


= 


h 




G2 


= edt 




h 





h 


= e 


'hi' - hex') 


G3 


= e^ 


dt 




h 


^ 


h 




G4 


= e2 


dx 




h 


= 


h 


= e2 


i 


G5 


= e2 


tdx 




h 


= 


h 


= e2 


(x — it) 


Ge 


= e2 \2tdt + 
= e^ \t^dt + 


xdx] 


k 


= 


h 


= e2 


i{x — it) 
i^{x - ti)^ 


G7 


txdx] 


h 




h 


1 

= 2 


Gs 


= e 


j^et^dt 
ht'dt 


+ {t+ ^et^x) dx] 


k 


= ieH-'x'' 


h 


= e 


ti - e {^t^ - ^t^xx + ht'x^)] 


G9 


= e 


+ (1 + et^x) dx] 


k 


= etx'^ 


I9 


= e 


i - e (te^ - t'^xi + ht^x'^)] 



(3.9) 

The only Noether symmetry which produces a useful integral is Gi. Had one treated 
only to 0(e), which is the implication of the one of the few texts giving a treatment 



of approximate symmetries [11|, one may have become enthusiastic about the integrals 
produced by G7 and Gs- Indeed they appear to be classic examples of a perturbation 
solution. However, we see from the O (e^) results that these symmetries are simply the 
result of truncation since the order of e increases with the order of the approximation. We 
remark that the results for this Lagrangian are not typical of Lagrangians of the general 
Emden-Fowler equation. It is well-known that the cases n = 2 and n = —3 differ from 



the general case |22|. Were we to use another value of n, the number of approximate 
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symmetries would differ and the symmetries, Gj and Gg, at the first order would not 
occur. We note that the Lie symmetry, G2, of this Emden-Fowler equation is not a 
Noether symmetry. 



4 Lie integrals 

The attraction of Noether's Theorem is that the procedure for obtaining a first inte- 
gral/invariant is trivial once the symmetry and gauge function have been determined, a 
process which is simultaneous. Of course there is a price to pay and that is the close 
relationship between the dependence of the coefficient functions on the derivatives and 



that of the integral/invariant |25|. One must look to generalized symmetries to be assured 
of finding the requisite number of invariants to know that a system is integrable. Apart 
from the practical problem of a more complicated calculation there is always the concern 
that one may have used the wrong Ansatz for the functional dependence on the derivatives 
in the coefficient functions. The Lie method is kinder in that respect, but the kindness is 
tempered with the additional burden of not completely friendly calculations. For multidi- 
mensional systems a combination of Noether's Theorem coupled with Lie symmetries has 
been found to be advantageous, particularly as, in the Lie approach, one can always insist 
that the integral be invariant under more than one symmetry Q. 

We noted at the conclusion of the previous section that the Lie symmetry, G2, listed 
in (|2.9|) is not a Noether symmetry and so does not have the easy route to an invariant 



which Noether's Theorem provides. We know from the existence of Ii in ( |3.9| ) that the 
Emden-Fowler equation (2^) is solvable as a consequence of Liouville's theorem |19]. 
However, can we look to G2 to provide an invariant which with Ii will give the solution? 
For this we use the Lie method outlined in § 1, ( |1.8| ). The invariants of G2 = td/dt—2xd/dx 
are found from the solution of the associated Lagrange's system 



dt dx dx 



(4.1) 



t —2x —3x 
and are 

u = xt^ and v = xt^ . (4.2) 
The second requirement of (|1.8[ ) is now expressed as 

dl{u{t,x),v{t,x)) ^ ^ 

dt ^ ' ' 

and leads to the second associated Lagrange's system, 
du dv 



V + 2u ev? + 3f 



(4.4) 



which is an Abel's equation of the second kind and, like most Abel's equations of either 
kind, is not responsive to known methods of integration. Consequently the existence of 
a second exact Lie point symmetry does not help in the search for a second integral. 
In terms of the conventional approach one is not surprised as the Lie Bracket of Gi 
and G2 is [Gi,G2] = Gi and we are effectively reducing the order by means of the non- 
normal subgroup, G2, which is well-known to remove Gi as a point symmetry of the 
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reduced equation. However, one may use an unconventional approach to attain our desired 
objective. In the reduction of order we know that the symmetry, Gi, becomes nonlocal. 
We know that the type of nonlocality is quite specific. The point symmetry becomes 
an exponential nonlocal symmetry and exponential nonlocal symmetries are known to be 
sufficient for the reduction of order [^. This is not precisely what we wish to do here 
since the second order ordinary differential equation has already been reduced to a first 
order ordinary differential equation and the further reduction to an algebraic equation, 
while perhaps of some mathematical interest and pleasure, does not really further for the 
solution process. What we wish to do is to use the nonlocal symmetry to transform the 
first order ordinary differential equation to a form such that the solution of the first order 
ordinary differential equation is reduced to an elementary procedure. The reversion of the 
process of reduction of order to a first order ordinary differential equation which can be 
solved is not normally a trivial task. Here this is not a concern as we already have one 
first integral and that coupled with the invariant we seek now is all that is required. (Of 
course we know that the solution of this Emden-Fowler equation is expressed in terms of 
elliptic functions and so we do not expect to be able to write x{t) in a simple form.) 
The symmetry, Gi, becomes 



exp 



du 



2u + v 



d d 
2u— + 2,v— 

ou ov 



(4.5) 



when we reduce the order of ( p.7[) to the first order ordinary differential equation (4.4). 
Were (4.5) a point symmetry we would seek a transformation which would make ( [4. 41) 
autonomous by converting ( |4.5| ) to a simple d/dX. However, ( ^.5] ) is not a point trans- 
formation and we cannot expect to be able to achieve that much. We use a modification 
of the standard method which seeks a transformation [^, p. 110] 

X = F{u,v), Y = G{u,v) (4.6) 

so that the functions F and G satisfy the linear first order partial differential equations 



XiF = 1 



and 



XiG = 0. 



(4.7) 



We simply replace the 1 in the first of (4/?) by the exponential part of the nonlocal 
symmetry (4.5), so that in the new coordinates Xi will still be an exponential nonlocal 
symmetry of the form 



X, 



exp 



du 



2u + v 



d 

dX' 



(4.8) 



Naturally the integral would be expressed in the new variables when the relationship is 
known. This stratagem means that the exponential term does not enter into any calcula- 
tions. The characteristic for the independent variable is easily computed to be 



w 



VU 2 . 



(4.9) 



The other characteristic (for the first of ( ^.T] )) requires the solution of 



(4.10) 



154 



P G L Leach, S Moyo, S Cotsakis and R L Lemmer 



which very easily gives 

C = F-hogu. (4.11) 
Under the change of variables 

X = vu~^, Y = logu^ (4-12) 



(4.4) takes the somewhat more friendly form 



dF X + 2exp[-y] 

dX 3X^-2e ^ ^ 

which is linear in the variable exp[y] and we immediately obtain the invariant 

J = e^sx^ - 2e'/' + 2 / (3^2 - 2e) dX. (4.14) 



That the expression for J is complicated is not surprising, but it is a precise form and 
maybe gives a different approach to the investigation of elliptic functions. 

The case of the Emden-Fowler equation has not given much joy for the use of ap- 
proximate symmetries to obtain approximate integrals/invariants as there are none which 
satisfy the criterion of not disappearing with increasing powers of the small parameter e. 
However, we can illustrate the way approximate symmetries can give truely approximate 
integrals/invariants with the simple harmonic oscillator. We should not be surprised at 
this for, if something does not work for the simple harmonic oscillator, it shall never work 
for anything. One regrets that there are some procedures which work only for the simple 
harmonic oscillator and that is the only model upon which trials are conducted! 



5 Discussion 

We have specified the approximate symmetries for the simple harmonic oscillator treated as 
a perturbation of the free particle, the Ermakov-Pinney equation treated as a perturbation 
of the equation for the simple harmonic oscillator and an autonomous Emden-Fowler 
equation of order two considered as a perturbation of the free particle. In the case of 
the simple harmonic oscillator the number of exact symmetries is eight. The number 
of approximate symmetries was found to be sixteen at 0(e), twenty-four at O (e^) and 
thirty-two at O (e^)- The homogeneity symmetry xd/dx occurs at all orders of e. The 
type of approximate symmetries that arise can be put into three categories: Those that 
are common to both perturbed and unperturbed equations; then the truely approximate 
symmetries the level of approximation of which increases with increasing powers of e. The 
third category consists of those symmetries in which every term is dependent on e. The 
Ermakov-Pinney equation has three exact symmetries instead of eight. When treated as a 
perturbation of the equation for the simple harmonic oscillator it has eleven approximate 
symmetries at 0(e), fourteen at O (e^) and seventeen at O(e^). The Ermakov-Pinney 
equation as a perturbation of the free particle has twelve approximate symmetries at 0(e), 
fourteen at O (e^) and sixteen at O (e^) . Another aspect which has been investigated is 
that of Noether symmetries and approximate first integrals (or invariants). 
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It is important to note that the knowledge of exact symmetries of a differential equa- 
tion enables us to reduce it and, if possible, completely solve it to obtain exact solutions. 
However, if one seeks for approximate rather than exact solutions, then approximate sym- 
metries are as useful as exact ones and can be determined perturbatively. In this paper 
we have emphasized more the study of approximate symmetries rather than approximate 
solutions. Finally we should note that in principle there exists the possibility of the intru- 
sion of secular terms into the expressions for the approximate symmetries and integrals. 
This is not the case with the equations and integrals treated here. As applications of 
approximate symmetries have so far been at the first order in e, whence the motivation of 
this paper which was to look beyond the first order in e, there is no evidence that secular 
terms and lack of convergence would present a problem. The purpose of this paper was to 
demonstrate the different types of approximate symmetry available, a result which is not 
accessible when approximate symmetries are treated to the first order only as has been 
the case in the earlier literature. 
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